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Abstract. The microphysics of relativistic collisionless sheared flows is investigated 
in a configuration consisting of a globally neutral, relativistic e“e+ beam streaming 
through a hollow plasma/dielectric channel. We show through multidimensional PIC 
simulations that this scenario excites the Mushroom instability (MI), a transverse 
shear instability on the electron-scale, when there is no overlap (no contact) between 
the e“e+ beam and the walls of the hollow plasma channel. The onset of the MI 
leads to the conversion of the beam’s kinetic energy into magnetic (and electric) 
field energy, effectively slowing down a globally neutral body in the absence of 
contact. The collisionless shear physics explored in this configuration may operate in 
astrophysical environments, particularly in highly relativistic and supersonic settings 
where macroscopic shear processes are stable. 
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1. Introduction 

Sheared flow conflgurations are pervasive in nature. Their study is of fundamental 
importance since they trigger instabilities and drive turbulence both in ionized and 
non-ionized settings. While the stability and evolution of sheared flow conflgurations 
has been thoroughly studied at the macroscopic (hydrodynamic/magnetohydrodynamic) 
level [B El El 51 El E!, their dynamics at the microscopic (plasma kinetic) level 
remain poorly understood. Recent theoretical and particle-in-cell (PIC) simulation 
studies have revealed that unmagnetized sheared flow configurations host a wealth 
of collisionless electron-scale instabilities m El El EDI [m 521 [IH]. These include the 
electron-scale Kelvin-Helmholtz instability (ESKHI), which develops parallel to the 
flow [a El El EDI EB E2], and the Mushroom instability (MI), which develops in 
the transverse direction to the flow [13]. Both the ESKHI and MI result from the 
coupling between high-frequency electromagnetic waves and electron plasma waves 
in the presence of a velocity shear. These microscopic instabilities can play an 
important role in ultra-relativistic and highly supersonic sheared flow scenarios, where 
the (magneto) hydrodynamic Kelvin-Helmholtz instability (KHI) is found to be stable 

HE]. 

Moreover, due to their electromagnetic character, both the ESKHI and MI can 
operate in the presence of a finite vacuum gap between shearing flows, i.e., in the absence 
of ’’contact” na. This effect is closely related to the development of electromagnetic 
instabilities in shearing, globally neutral, polarizable metallic/dielectric slabs separated 
by a nanometer-scale gap, which results in an effective non-contact friction force between 
the slabs OEH]. This form of non-contact friction, associated with the development 
of electromagnetic instabilities, is the classical analogue to the quantum friction effect 
proposed by Pendry [E|. While the effect of non-contact friction has been considered 
in the sub-relativistic regime, connected to the development of electromagnetic modes 
parallel to the flow (ESKHI), the role of electromagnetic modes transverse to the flow 
has been overlooked. In this work, we show that the transverse modes are dominant 
in relativistic settings, and are responsible for the non-contact friction force in the 
relativistic regime. 

Reproducing in the laboratory the collisionless conditions that characterize 
most astrophysical systems is highly challenging. Sheared flow experiments have 
been performed in the collisional regime mi EE], revealing the development of 
the hydrodynamic KHI under high-energy-density conditions. More recently, the 
experimental observation of sheared flow driven turbulence in the collisionless regime 
was reported [T9|, but clear signatures of the ESKHI/MI were not observed. We believe 
that the ESKHI/MI processes may be better observed in configurations where bulk 
overlap between the flows is minimized, inhibiting streaming instabilities that could 
mask their effects. Relativistic, low divergence and globally neutral plasma flows may 
be achieved by superimposing relativistic electron (e“) and positron (e"*") bunches from 
linear accelerators [20]. Alternatively, the generation of energetic, globally neutral 
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beams has also been achieved resulting from QED pair cascades triggered by the collision 
of a laser-wakefield-accelerated electron bunch with a solid target 1211 ^ The production 
of such exotic beams will enable the laboratory exploration of collisionless microphysical 
processes in the relativistic regime, which are directly relevant to many astrophysical 
systems. In particular, it has been suggested that the interaction of such beams with 
a bulk stationary plasma would allow to probe the collisionless current filamentation 
instability (CFI) in the relativistic regime [201 El]- Here, we investigate a configuration 
where an beam streams through a stationary hollow plasma channel, or hollow 

dielectric, in order to probe the microphsyics of relativistic collisionless sheared flows. 

In this work, we investigate the relativistic collisionless shear dynamics associated 
with a relativistic, globally neutral beam streaming through a hollow plasma 

channel. The e“e''‘ beam is considered to be propagating with a relativistic Lorentz 
factor 7 o, and to have a gaussian density profile with a± = 2a\\ = 2clujpe (T and || 
are perpendicular and parallel to the propagation direction, respectively), and a peak 
number density ni, (measured in the laboratory frame); c denotes the speed of light in 
vacuum, and ojpe = y^^e^nbjmfef) where \/2ujpf, is the plasma frequency of the beam (the 
factor arises from the positron contribution). These beam parameters have been 
considered in the context of exploring the CFI in the relativistic regime |20|, and can 
be made available at SLAG National Accelerator Laboratory with 70 ~ (9(10^). The 
hollow plasma channel has an internal radius R where the plasma density is assumed to 
be 0, and the surrounding plasma density is Up = n^,. This channel could also be made 
of a hollow solid-state dielectric. The impact of the dielectric properties of the hollow 
channel on the development of the shear dynamics will be explored elsewhere. In this 
work, we assume the channel walls to have the dielectric response of a cold collisionless 
plasma with the same density as the beam. 

We first explore the changes in the dynamics of the propagation of the relativistic 
beam in channels with different radii R. In particular, we aim to clearly 
distinguish the dynamics in narrow hollow channels (i^/cr^ ~ 1 ), where significant 
overlap between the beam and the walls of the hollow channel occurs, triggering bulk 
streaming instabilities like the current filamentation instability (CFI) [201 121] , from the 
dynamics in wide hollow channels {R/a±^ > 3), where bulk overlap is negligible and 
only collisionless sheared flow instabilities may develop. We thus begin by reviewing the 
growth rates of the relevant instabilities in order to understand their relative importance 
under different conditions. 

2. Bulk streaming instabilities versus shear instabilities 

The the development of the CFI triggered by a relativistic beam propagating in 

a stationary bulk plasma has been explored by [201 ED- The CFI growth rate in the 
relativistic limit scales as rcFi/<^pe ~ H is expected that if significant overlap 

between the beam and the walls of plasma channel exists, the CFI will develop with 
Tcfi calculated using the average overlapping beam density. 



Slow down of a globally neutral relativistic e e+ beam shearing the vacuum 4 

We have shown elsewhere [aim US! that unmagnetized collisionless sheared flows 
can trigger both the ESKHI (along the sheared flow) and the MI (transverse to the 
sheared flow). Although the f ini te geometry of the beam will certainly be important, we 
can estimate the growth rates of both MI and ESKHI modes based on a slab geometry 
for the sheared flow. Assuming a semi-infinite relativistic plasma slab propagating with 
7o, shearing with a second semi-infinite stationary plasma slab, the respective growth 
rates scale as EESKHi/^pe - 7o 0 E] and rMi/<^pe ~ 7 o [IS]- It is therefore 
expected that under highly relativistic conditions, relevant to the e“e+ beam streaming 
in a hollow channel scenario, the transverse MI will be dominant due to its weaker 
dependence on 70 . In addition, the length of the e“e''‘ beam (~ a\\ = cjiVpf) is much 
shorter than the wavelength of the most unstable mode of the ESKHI (~ so 

the ESKHI should not develop at all. Eurthermore, due to their electromagnetic nature, 
these microscopic shear instabilities can operate in the absence of ’’contact”, i.e. when 
there is a finite gap of vacuum between the sheared flows ra. which can be directly 
explored in the hollow plasma channel configuration. Using the appropriate boundary 
conditions to describe the shear interaction between an ultra relativistic and stationary 
plasma slabs, separated by a finite gap Lg, one finds that the growth rate of the MI 
scales as l/A^. 

3. PIC simulations of a relativistic e“e''‘ beam streaming through a hollow 
plasma channel 

We analyze the effect of different diameters of the hollow plasma channel on the 
dynamics of relativistic beam. We have performed 2 D PIC simulations of the 

transverse plane to the propagation direction using the PIC code OSIRIS [221 1^ . 
Although the 2 D simulations assume both an infinitely long beam and hollow plasma 
channel, these simulations capture the transverse electron-scale shear dynamics and the 
interplay with the CPI in the case of overlap. We simulate a domain with dimensions 
L, X Ly = 26^ {cjojpeY, resolved with 256^ cells; 36 particles per cell per species were 
used. The beam is positioned at the centre of the box and propagates out of 

the plane (- 1 -;^ direction). The walls of the cylindrical plasma channel are placed at a 
distance R from the centre of the beam. The plasma composing the walls of the hollow 
plasma channel has the same electron density as the beam, and contains heavy ions 
(assumed immobile in the simulations). 

3.1. Strong overlap 

Pigure illustrates the case of strong overlap between an e“e+ beam (with 70 = 10 ^) 
and the walls of a narrow hollow plasma channel with R = a^ = ^cjujpf.. The evolution 
of the positron density of the beam and the electron density of the walls of the 

plasma channel is shown in Pigure[T](a), and the evolution of the self-generated magnetic 
field is shown in Pigure[^ (b). The strong overlap between the beam and the walls of the 
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channel is clearly observed in Figure The initial condition of current neutrality 

is verified in the magnetic field structure of Figure (bl), which shows small amplitude 
and small-scale magnetic fluctuations due to thermal noise of the beam. At Upet = 80 
the development of magnetic field structures associated with the CFI becomes clearly 
visible (inset of Figure (b2)); these structures occur precisely in the region of overlap 
between the beam and the walls of the plasma channel, creating the bulk counter¬ 
streaming conditions required to trigger the CFI [20]. The growth rate of the CFI, 
measured from the evolution of the total magnetic field energy in the simulation domain, 
matches the theoretical growth rate when taking into account the density of the beam 
in the overlapping region. The nonlinear development of the CFI ultimately leads to the 
break-up of the beam into multiple current filaments as shown in Figures (a3) 

and (b3). 

3.2. Weak overlap 

For channels with R > 4a±, the overlapping beam density with the walls of the 
plasma channel becomes very weak, inhibiting the development of the CFI. A case 
of weak overlap is illustrated in Figure where the channel radius has R/cr± = 4. All 
other parameters are kept fixed relative to the case of Figure [T] Again, the system is 
initially current and charge neutral, and no background fields are present. Only small 
electromagnetic fluctuations are present in the early stage of the simulation, associated 
with the finite emittance of the relativistic beam. At iVpet = 600, a self-generated 
magnetic dipole structure emerges as shown in Figure (b2). This is a signature of 
the transverse electron-scale shear instability, the MI. For small a±, as used in these 
simulations, the dipole mode is found to be the fastest growing mode. Higher order 
modes have been observed in simulations with larger values of a±, while fixing R/a± = 4. 
A predictive analytic theory describing this configuration is still lacking, and will be 
developed in a future publication. The continued growth of the dipole magnetic field 
structure due to the MI, eventually leads to the breakup of the beam (Figures 

(a3) and (b3)). If the beam escapes the channel before breakup, the dipolar magnetic 
field structure can be preserved and sustained by the beam on a time-scale 1 / (cOpg/ . 
This configuration can be used to imprint a tailored field structure in a globally neutral 
relativistic beam. 

The growth rate of the total self-generated magnetic field energy (measured from 
2D simulations) for different hollow channel radii (R) is shown in Figure]^ (a). The cases 
with R/o'± < 3 are characterized by strong overlap between the beam and the walls of 
the hollow channel. The growth rates are found to agree with the theoretical prediction 
for the CFI (dashed line), taking into account the beam density in the overlapping 
region. At R/cr± = 3, magnetic field generation activity occurs both in the overlapping 
region, due to the CFI, and at the beam position, due to the MI, marking the transition 
to the collisionless shear-dominated regime. For R/o'± > 4, the measured growth rates 
are higher than predicted by the CFI curve, which decays abruptly due to the low beam 
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density at overlap. Here, the dynamics of the system is dominated by the collisionless 
shear dynamics of the MI, whose growth rate is found to decay as {R/cr±)~^/^. 

Figure]^ (b) shows the dependence of the MI growth rate on the relativistic Lorentz 
factor ( 70 ), for a channel with R/a± = 4, taken from 2 D and 3D PIC simulations. The 
growth rates measured from 3D simulations (blue points) are ~ 5 times lower than 
those measured from 2 D simulations (red points). This numerical factor arises from the 
spatio-temporal development of the MI in the longitudinally finite (cr|| = Ic/iVpe) beam 
in 3D. Nevertheless, the MI growth rate is found to decay in both cases with 7 q ' , 
which is consistent with the theoretical model for the slab geometry da. 

The fully self-consistent structure of the Ml-generated magnetic fields from a 3D 
PIC simulation is illustrated in Figure The simulation domain has dimensions 
L^xLyXLz = 26^{c/u)peY, resolved with 256^ cells. The red/orange isosurfaces represent 
the positron density of an ultra relativistic e“e+ beam with 70 = 10 ^, propagating in the 
-|-x-direction inside a hollow plasma channel (grey isosurfaces) with R/o'± = 4. Half of 
the plasma channel wall has been omitted to observe the growing magnetic field structure 
inside. The Ml-generated dipolar magnetic field structure, previously identified in 2 D 
simulations, is clearly observed at the beam position in the 3D simulation (Figure]^ 

(b) ); the field lines are bowed towards the back of the beam in 3 D (Figure]^ (a)). The 
growing magnetic field is sustained by a current structure associated with the charge 
separation between electrons and positrons of the beam. This charge separation also 
induces a dipolar electric field structure, which triggers the trailing wake-field structure 
behind the beam (Figure]^ (a)). For an beam density of ni, = 10^^ cm“^, the snap 
shots of Figure [^correspond to a propagation length of ~ 50 cm. 

The generation of electric and magnetic field energy due to the MI is accompanied 
by the self-consistent loss in kinetic energy of the beam. This is illustrated in Figure [^ 

(c) , which shows the exponential growth of electric and magnetic field energy, and the 
corresponding negative variation of the beam’s kinetic energy. The development of 
the MI therefore acts as an effective non-contact friction force, which is proportional 
to the MI growth rate, dp\\/dt = — 2rMi P\\ = Dn.c. {p\\ is the average longitudinal 
momentum of the beam). The effective non-contact friction force (D^ c ) associated with 
the development of electromagnetic instabilities between shearing dielectrics separated 
by a vacuum gap has been previously explored in m [T5] . where the longitudinal 
electromagnetic modes played the most important role at subrelativistic velocity shears. 
Here, we have shown that the transverse MI modes are responsible for the effective 
non-contact friction effect in the relativistic regime. 

3.3. Non ideal e“e+ beam and off-axis propagation 

If the e“e+ beam is not perfectly charge neutral (ne+ — n^- = 5n 0) it will defocus 
on a time-scale Wpe'Pjefocus 3 o^‘^inb/dnyR, where is min(ne+,Ue-). Fortunately, for 
ultra relativistic beams, this defocusing time is much larger than the growth time of 
the MI LVpeTui ~ <^pe/rMi oc 7 ^^, allowing for the development of the MI before the 
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beam defocuses. We have observed the persistent development of the MI in beams 
with 5n = O.lrift in 3D PIC simulations, where the Ml-generated fields were found to 
be superposed over the self-consistent fields associated with the global cnrrent of the 
non-nentral beam. 

We have also simulated the effect of the e“e''‘ beam propagating off-axis, at a finite 
distance vq from the axis of the hollow channel. We find for offsets oiro/R < 1/4, while 
still preserving the weak overlap condition (R — ro)/a± > 4, the continual development 
of the MI with only small variations in the growth rate < 10% relative to the ideal 
on-axis case. The emerging dipolar field structure becomes slightly asymmetric under 
these conditions. 

4. Conclusion 

In conclnsion, we have shown that the transverse microphysics of relativistic collisionless 
sheared flows can be investigated in a scenario where a relativistic, globally neutral 
beam propagates inside a (stationary) hollow plasma channel. Resorting to 2D PIC 
simulations, we have stndied the propagation of a relativistic e“e+ beam in different size 
hollow plasma channels. We have found that for channel radii with R/a± < 3 strong 
overlap between the channel and beam occurs, exciting the CFI in the overlapping 
region, whereas for channel radii with R/a±> 4: the overlap is weak, inhibiting the CFI 
and allowing for the clear development of the relativistic collisionless shear dynamics 
of the MI. For the range of e“e+ beam and hollow channel parameters explored in this 
work, the MI growth rate was found to scale as PMi/<^pe oc Interestingly, 

despite the microscopic natnre of this shear instability, we find the persistent onset of 
the MI even when the gap between the beam and the walls of the plasma channel are 
R/<7± S> 1. These scalings also support the development of the MI when the e“e+ 
beam is not perfectly neutral, since the MI growth time is shorter than the electrostatic 
defocusing time. The full spatio-temporal evolution of the MI was observed via 3D PIC 
simulations, having revealed the self-consistent strncture of the Ml-generated magnetic 
field. The self-consistent slow down of the relativistic e“e''‘ beam, due to the conversion 
of the beam’s kinetic energy into electric/magnetic field energy through the MI was also 
verified. The development of the MI therefore leads to an effective non-contact friction 
force which is proportional to Fmi- These are the first ab initio 3D PIC simnlations 
capturing the self-consistent slow down of a globally neutral body in the absence of 
contact [m [T51 [TB] . 
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Figure 1. Evolution of (a) number density and (b) self-generated magnetic field 
in the interaction between a relativistic e“e+ beam (70 = 100 , propagating out of 
the plane) and a narrow hollow plasma channel (at rest) with R = cf±, illustrating 
the development of the bulk CFI (inset of b 2 ) due to the strong overlap between the 
beam and the plasma channel. Red and grey scales correspond to the positron number 
density of the beam and the electron density of the hollow plasma channel, respectively. 
Evolution of the interaction is shown at (1) ujpet = 20 (initial setting), (2) ujpet = 80 
(linear regime) and (3) Upet = 140 (nonlinear regime/beam breakup). 
































Slow down of a globally neutral relativistic e e+ beam shearing the vacuum 


10 


ne+ beam [Hq] cylinder [Hq] 
0 1 2 0 1 2 



I 10 


3 

:6 I 

'h 

- 4 ^ 
pq" 
^2 m 


LJO 

■ 10 


6 B 
o 

4 ^ 
pq" 
2 m 


Uo 

25 


I 


20 ^ 
3 


15 

10 

5 

LJO 


Figure 2. Evolution of (a) number density and (b) self-generated magnetic field in the 
interaction between a relativistic e“e+ beam (70 = 10 ^, propagating out of the plane) 
and a wide hollow plasma channel (at rest) with R = (negligible overlap between 
beam and plasma channel), illustrating the development of the transverse relativistic 
shear instability (MI). Evolution of the interaction is shown at (1) ujpet = 50 (initial 
setting), (2) ujpet = 600 (linear regime) and (3) Up^t = 850 (nonlinear regime/beam 
breakup). 
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Figure 3. Left: growth rate dependence of the self-generated magnetic field on the 
inner radius R of the hollow plasma channel, showing transition from CFI to ML Here, 
the e“e+ beam has 70 = 10^ and a fixed transverse width of = 2clijjpe. Red points 
represent results of 2 D simulations. Dashed curve represents theoretical growth rate 
of the CFI considering the beam density at overlap. Right: Growth rate of the MI- 
generated magnetic field versus e“e+ beam Lorentz factor 70 . Here, the inner radius 
of the hollow channel is set to R = 4(7_l (weak overlap/no CFI), with a± = 2clujpe. 
Red and blue points represent results of 2 D and 3D PIC simulations respectively. 



Figure 4. 3D PIC simulation of the development of the MI for 70 = 10^, 
(Jx = 2.(7\\ = 2cju^Q and R = Lax- The e“e+ beam density is represented by the red 
isosurfaces, and the hollow plasma channel is represented by the grey isosurfaces (half 
of the cylindrical structure has been omitted to more clearly illustrate the dynamics of 
the e“e+ beam inside). The e“e+ beam propagates in the positive x direction. The red 
lines represent the structure of the Ml-generated magnetic field. Figures (a) and (b) 
are side and front views of the simulation domain, respectively, taken at UrpR. = 36500. 
Figure (c) shows the temporal evolution of the electric (orange curve) and magnetic 
(blue curve) field energy and the negative variation of the total kinetic energy of the 
beam (red curve), all normalized to the initial kinetic energy of the beam. The slope 
of the dashed line represents the growth rate of the MI. 

















